A new broken [/(l)-symmetry in extreme type-II superconductors 
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A phase transition within the molten phase of the Abrikosov vortex system without disorder 
in extreme type-II superconductors is found via large-scale Monte-Carlo simulations. It involves 
breaking a !7(l)-symmetry, and has a zero-field counterpart, unlike vortex lattice melting. Its 
hallmark is the loss of number-conservation of connected vortex paths threading the entire system 
in any direction, driving the vortex line tension to zero. This tension plays the role of a generalized 
"stiffness" of the vortex liquid, and serves as a probe of the loss of order at the transition, where a 
weak specific heat anomaly is found. 
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The key role of topological excitations in all phase- 
transitions, whose proliferation is accompanied by the 
loss of a generalized "stiffnessL lias been emphasized ku. 
Andersoru. Recent numericaltTEl and analytical workO'El 
reveals that in extreme type-II superconductors, topolog- 
ical excitations in the form of vortex loops are essential 
to the physics of the vortex system in low magnetic fields. 

It will be shown in this paper that the proliferation of 
large vortex loops induces a phase transition within the 
molten phase of the Abrikosov vortex system without dis- 
order. The phase-transition thus appears in a region of 
the phase-diagram where one normally would not expect 
critical behavior simply based on a study of the local 
superconducting order parameter ^ entering Ginzburg- 
Landau theory. It is, however, a general property of 
extreme type-II superconductors, distinct from the first 
order vortex lattice melting transition in that it has a 
zero-field counterpart, i.e. the transition from a super- 
conductor ptai-a-iiprrnal metal mediated by a vortex-loop 

"biowout"&mM 

The model for extreme type-II superconductors con- 
sidered in this paper, is the uniformly frustrated 3 di- 
mensional anisotropic XY (3-DAY) modelEJ, defined by 
the Hamiltonian 



H({8(r)}) = - J » C0S I V ^(r) - 



r,fj,—x,y.z 



(i) 



where 9 is the local phase angle of the superconduct- 
ing complex order parameter and V is a lattice deriva- 
tive. The coupling energy along the /i-axis, J^, is de- 
fined by J x = J y = (<S> 2 Q d)/(16Tr 3 \ 2 ab ) = J ± , and J z = 
(<I>Q^ fc )/(167r 3 A^d), where A a & and A c are the magnetic 
penetration lengths set up by screening currents in the 
CuO plane and along the crystal c-axis, respectively. $o 
is the flux quantum, £ a b is the superconducting coher- 
ence length within the CuO-planes, and d is the distance 
between two CuO-layers in adjacent unit cells. We may 
use the lattice spacing as a measure of In Eq. |l|, 
is related to the non-fluctuating external magnetic vec- 
tor potential A vp by A^r) = 27r/$ J r r+e " 
where e M is the unit vector along the /x-axis 



We consider systems of size L x x L y x L z . To per- 
form simulations and finite size scaling of systems with 
very low filling fractions, we define an "extended" Lan- 
dau gauge 
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, m are positive integers satisfying n x n 



dr' ■ A vp (r'), 



where n x , n. 

L y , and m y m = L x . Hence, the filling fraction / is given 
by / = n m [n x - m y ]/L x L y . 

In this paper, we consider the filling fractions / in the 
range 1// G (20, .., 1560). These filling fractions are so 
low that spurious pinning to the numerical lattice is elim- 
inated, i.e. any spurious transverse Meissner- effect has 
vanished, well below all temperatures of interest in this 
papeM-R Using / = B& b /<f> and U = 151, 1// = 1560 
corresponds to a uniform induction B ~ 0.56T. Using 
Kit > Su as a condition for uniform B and X a b — 1500A, 
we find that this assumption is valid for inductions as 
low as 0.1T, corresponding to a filling fraction of or- 
der 10 -4 . The uniform induction B is taken along the 
crystal c-axis. The anisotropy parameter is given by 
T = V ' J_l / ' J z — \ c d/ \ a b$,ab- For each filling fraction, 
we carry out simulations on three systems with different 
sizes (~ 40 3 ,~ 80 3 ,~ 150 3 ) to study finite -size effects. 
For the specific heat, we also consider 360 3 for / = 1/90. 

In addition to the specific heat CoQ, we calculate three 
quantitities in this paper. The first two are standard, 
while the third is unusual and probes a subtle U(l)- 
symmetry of the system. 

i) Helicity modulus. To probe the global phase co- 
herence we consider the helicity moduli, T M , defined as 
the second derivative Qf-ithc free energy with respect to 
an applied phase twistOH When T M is zero, the resis- 
tance along /Lt-direction is finite, and any applied current 
along /i-direction will dissipate energy. Note that T M is a 
global quantity, so even when T M = the system can still 
maintain local phase coherence and exhibit diamagnetic 
fluctuations. T z will vanish at the temperature T z , coin- 
ciding with the zero-field superconducting transition T c 
when / — 0, and will coincide with the melting temper- 
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ature of the Abrikosov vortex lattice (AVL) when / > 0. 

ii) Structure factor. To probe the Abrikosov vortex 
lattice (AVL) melting, wexonsider the structure function 
for q z vortex segmentsEm. A vortex configuration can 
be obtained from a phase-configuration from the coun- 
terclockwise line integral of the gauge-invariant phase- 
differences around any plaquette of the numerical lat- 
tice. It must satisfy the condition for conservation of 
vorticity, J2 e [V v 6(r) - A v (t)] = 2^(r) - /„). Here, 
C is the closed path traced out by the links surrounding 
an elementary plaquette, and v represents the Cartesian 
components of the links comprising the closed path C. 
Furthermore, ^(r) = 0, ±1 represents a vortex segment 
penetrating the plaquette enclosed by C. If we moni- 
tor the structure function 5(Q, k z = 0) of one chosen 
Bragg peak as a function of temperature, we expect that 
5(Q, k z = 0) has a discontinuous drop to a very small 
value at the melting temperature. The structure func- 
tion will vanish discontinuously at the melting tempera- 
ture T m (B) of the AVL. 

Hi) Vortex-path probability Ol- We define Oi as the 
probability of finding a directed vortex path threading 
the entire system transverse to the induction B, without 
using the PBC along the field direction. It is obtained 
by computing the number Ny of times we find at least 
one such path threading the system in any direction _L B 
in Np different phase-configurations, normalized by Np, 
i.e. Ol = Ny/Np. Note that this differs from algorithms 
employed in previous publicationsEjO. 

Ol = means that there is no connected path of vor- 
tex segments that threads the entire system in the trans- 
verse direction, without using PBC along the field direc- 
tion several times. Now, let N" (a € [x, y, z\) denote the 
areal density of connected vortex paths threading the sys- 
tem in any direction, including the direction parallel to 
the induction. It is clear that in the AVL phase Ol = 0, 
and N[ = B/$ , while N£ = N V L — 0. Thus, A£ is a 
conserved quantity at fixed induction B. On the other 
hand, Ol = 1 implies that N^' y > 0, and the total num- 
ber of vortex paths threading the system in any direction 
scales with system size, but undergoes thermal fluctua- 
tions. Therefore, N£ is no longer a conserved quantity. 
The change in Ol takes place at the temperature Tl(B), 
which coincides with T c when / = 0. 

Number conservation uniquely identifies a U(l)- 
symmetry, and hence the low-temperature phase of the 
vortex-system (the dual of the phase-representation of the 
superconductor) exhibits explicit [/(l)-symmetry, since 
Ol = 0. At high temperatures Ol = 1, is not con- 
served, and the [/(l)-symmetry is broken. Under no cir- 
cumstance can a [/(l)-symmetric phase be analytically 
continued to a [/(l)-nonsymmetric one. The change in 
Ol from to 1 therefore signals a geometric transition 
taking place withing the vortex-liquid, indicative of the 
breakdown of vortex- line liquid picture of the molten 
phase of the Abrikosov vortex lattice. 

We first discuss the zero- field simulation results. In 



Fig. we plot the specific heat C, the helicity moduli 
T X ,T z , and Ol as functions of temperature for a system 
with Size = 140 3 and anisotropy r = 7. C has a near- 
logarithmic singularity at ksTc/ J±_ = 1.12, where also 
and T z drop to zero. This second order phase tran- 
sition is the superconducting-normal state phase transi- 
tion. Note that both T z and T x vanish at T c , albeit with 
different amplitudes due to the anisotropy of the model. 
Since there is one, and only one, transition in zero field, 
this is a check that the system sizes used m the simula- 
tions are adequate for the given anisotropic!. 
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FIG. |^. Specific heat, vortex-path probability Ol, and 
helicity moduli along ^-direction (T x ) and along ^-direction 
(T z ) as functions of temperature for the zero field case. V = 7, 
and system size is 140 3 . For T < T c , O l = 0. For T > T c a 
vortex loop blow out has taken place, i.e. Ol = 1. 

For T < T c , it is clear that L = and A£ is 
conserved and equal to zero. If for T > T c , we have 
O l = 1, then A£ > 0. This feature of O l is pre- 
cisely what we find in our simulations, see Fig. [|. This 
means that the system undergoes a 3DXY-transition, 
with a low-temperature U (l)-symmetric state and a high- 
temperature {/(l)-nonsymmetric state. Note that in 
terms of the ordinary superconducting orderparameter 
^ in Ginzburg-Landau theory in zero field, the situa- 
tion is the opposite as far as t/(l)-symmetry and sym- 
metry breaking is concerned: The more familiar phase- 
representation of the superconductor is related to the 
vortex-representation via a duality transformation which 
interchanges low and high temperaturestJ. 

Next, we consider finite magnetic fields. In Fig. top 
panel, we plot the specific heat fL the structure factor S, 
Ol and the helicity moduli T x tB and T z as functions of 
temperature for a system with filling fraction / = 1/90, 
Size = 72 x 80 x 80 and T = 7. We see from the top 
panel of Fig. that at fcsT!,„/J^ = 0.49 the structure 
factor 5* drops discontinuously from 0.2 to 0, indicating 
a first order AVL melting transition at T m . At T = T z> 
T z drops discontinuously from 0.55 to 0. We see that in 
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our simulations T z — T m , indicating that the AVL melts 
directly into an incoherent vortex liquid. Thus, in the 
vortex liquid ptiaa c. there is no global phase coherence 
in any directiorJaa'tLKiis conclusion also holds for the 
isotropic case T = ltLZHIa. 
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FIG. . Top panel: Specific heat C, structure factor S, 
vortex-path probability Ol and helicity moduli Y^, T z as 
functions of temperature for a system with a filling fraction 
/ = 1/90 and Size = 72 x 80 x 80. The melting tempera- 
ture is given by fcsT m /Jx = 0.49. At Tl{B), we also show 
C for a system of size 360 3 , illustrated by solid squares. The 
curve is shifted up by the amount 0.2 ks from the curve for i 
72 x 80 x 80, for clarity. Note the increase in the weak anomaly 
at Tl(B), which is the hallmark of a thermodynamic phase- 
transition. 

Bottom panel: Helicity modulus along the field direction T z 
and Oz, as functions of temperature for three systems with 
/ = 1/90 and sizes 36 x 40 x 40, 72 x 80 x 80, 144 x 160 x 160. 



In the top panel of Fig. , the specific heat shows a spike 
precisely at T m , indicating a discontinuous jump in the 
internal energy, and thus proving the existence of a first 
order phase AVL melting transition with a latent heat at 
T m . We now focus on the temperature range where Ol 
changes from to 1. For the filling fraction / = 1/90, Ol 
rises from to 1 , and reaches 1 at T = Tl well separated 
from T m and the crossover temperature Tb C 2, see top 
panel of Fig. . vortex loop blowout takes place. For in- 



creasing system size the width of the transition decreases 
rapidly, see the lower panel of Fig. . Note also how the 
low-temperature tail of Ol is suppressed when increas- 
ing the system size, while the temperature Tl{B) where 
Ol reaches the value 1 stays almost fixed, moving slightly 
down with increasing system size. Thus, in the thermo- 
dynamic limit, there exists a well defined temperature 
Tl where Ol rises sharply from to 1. This is found 
for all the filling fractions considered. For f = 1/90, 
we have also performed simulations for systems of size 
360 3 , to bring out the scaling of the additional weak spe- 
cific heat anomaly at Tl{B), see the top panel of Fig. 
. It shows conclusively that Tl{B) is a thermodynamic 
phase-transition. 

The position of the specific heat anomaly at Tl{B) for 
the system of size 360 3 appears at a slightly lower tem- 
perature than those where Ol starts to dip down from 
1 in smaller systems. This is because the O^-curves be- 
come sharper as L increases, and the position of the C- 
anomaly is at a position which may be estimated to be 
the limiting temperature at which Ol changes abruptly 
from to 1. Note that in the specific heat, the anomalies 
both at T m (B) and Tl(B), both equal in magnitude, are 
well beyond the noise-level in the simulations. For the 
other parts of the C-curve, the uncertainties arc of order 
the symbol size. 

The fact that Ol = 1 above Tl(B), indicates that the 
long-wavelength vortex-line tension e(T) vanishes. We 
may infer this from the change in Ol- 

Furthermore, we have checked how Tl{B) varies with 
a change of aspect ratio L x /L z of the system. For / = 
1/380, r = 7, (not shown) we have considered systems of 
size L x x L y x L z /a, with L x = L y = L z = 40,80, 120, 
varying a £ [1.00, 1.25, 1.50, 1.75, 2.00]. Within a vortex- 
line liquid picture, we should find Tl oc oS3. Instead, we 
find a change of less than 5% on increasing a from 1 to 2, 
again difficult to explain within a vortex-line liquid pic- 
ture of the molten phase above Tl(B). This is support- 
ive of a change in the connectivity of the vortex-tangle 
at Tl(B), in a direction perpendicular to the magnetic 
field. In a vortex liquid with non-zero s(T), the vortex 
system is only connected across the system along the field 
direction. 

The long-wavelength vortex-line tension may thus 
serve as a probe of the loss of order at the transition, 
playing the role of a generalized "stiffness" EJ characteriz- 
ing the two vortex liquid phases above and below Tl(B). 
More precisely, the two vortex-liquid regimes are charac- 
terized by breaking a [/(l)-symmetry of the dual the- 
ory of the Ginzburg-Landau theory, on crossing the line 
Tl(B) from below. 

We have obtained C, T 2 , S, and Ol for a wide range 
of filling fractions, see Fig. . The resulting (_B,T)-phase 
diagram for the uniformly frustrated 3DXY model with 
the anisotropy parameter r = 7 and B || c, is shown in 
Fig. . In zero field, a 3-DXF-transition separates the su- 
perconducting phase from the normal phase. In a finite 
magnetic field, we have three distinct phases: the AVL 
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for T < T m (B), the vortex liquid with line tension for 
T m (B) < T < T L (B), and the vortex liquid without line 
tension for T > Tl(B). Tl(B) is a critical line, the finite 
field counter part of the zero field vortex loop blowout at 
T = T c . In the low field regime, / < 1/600, the AVL 
melts directly into a vortex liquid with no linetension. 

The dimensionless criterion determining the low-fi 
regime of the melting line T m (B), is that fT 2 << ll 
For our case, the merging of T m (B) and Tl(B) oc- 
curs at /r 2 w 1/13, which is well within the low- field 
regime. Note that at more elevated fields, fT 2 > 1/4, 
the vortex lattice melting line is well described by both 
the XF-model and the 2D boson-analogy of the vortex- 
system, as-*ecently nicely elaborated on by Koshelev and 
NordborgEl However, from their Fig. 1, we note that 
deviations from universal properties of linelike melting of 
the vortex lattice starts to appear for fT 2 < 1/5, consis- 
tent with the above picture. 

T m (B) meets Tl{B) at a tricritical point, see Fig. . 
The continuation of the melting line below the tricriti- 
cal point is first order. We have observed a substan- 
tial increase in the ^-function spike at the melting tran- 
sition below the tricritical point, since the entropy in 
the transition— at Tl(B) now contributes to a first or- 
der transition!!!!. Note also that well above the tricritical 
point, the position of T m (B) is adequately described by 
a Lindemann-criterion applied to the -London-model in- 
cluding only field-induced vortex lined 2 .]. It is not clear, 
however, that such a model would give T 2 = in the 
entire molten phase. 
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FIG. . The intrinsic (B,T)-phase diagram of an extreme 
type-II superconductor, obtained from the uniformly frus- 
trated 3DXY model with r = 7 and B \\ c. k B T c /J ± = 1.12. 
The filling fractions / that are used are given by 1/f — 
90, 132, 290, 380, 506, 650, 1560. The lines are guides to the 
eye. The H C 2-\ine is found from the maximum value of the 
broadened peak of the specific heat in a finite magnetic field. 



In summary, we have used large-scale Monte- Carlo 
simulations to analyze the characteristics of vortex-paths 
in a model of extreme type-II superconductors in the ab- 
sence of disorder. An abrupt change in this characteris- 
tics reveals a phase-transition involving the breaking of a 
{/(l)-symmetry, which is suggested to exist in zero as well 
as finite magnetic field. It is therefore distinct from the 
vortex lattice melting transition. This symmetry break- 
ing is a consequence of a thermally driven vortex-loop 
"blowout" Ej. It results in a nonconserved number of ther- 
mally fluctuating connected vortex paths threading the 
entire system in any direction, including the direction 
perpendicular to a transverse magnetic field. In other 
words, the connectivity of the vortex-tangle changed, 
and the insensitivity of 'Tl(B) strongly suggests that this 
is a feature that survives in the thermodynamic limit. 

In principle it is also possible for a J7(l)-symmetry 
breaking to be first-order, although in zero magnetic field 
this only happens for small values of kEj. In finite field 
this could change. More work is clearly needed to in- 
vestigate in detail the universality class of the proposed 
transition at Tl(B). 
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